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ABSTRACT 

The wing leading edge and nose cone of the Space Shuttle are fabricated from a reinforced carbon/carbon material 
(RCC). The material attains its oxidation resistance from a diffusion coating of SiC and a glass sealant. During re-entry, the 
RCC material is subjected to an oxidizing high temperature environment, which leads to degradation via several 
mechanisms. These mechanisms include oxidation to form a silica scale, reaction of the Si0 2 with the SiC to evolve 
gaseous products, viscous flow of the glass, and vaporization of the glass. Each of these is discussed in detail. 

Following extended service and many missions, the leading-edge wing surfaces have exhibited small pinholes. A 
chloridation/oxidation mechanism is proposed to arise from the NaCl deposited on the wings from the sea-salt laden air in 
Florida. This involves a local chloridation reaction of the SiC and subsequent re-oxidation at the external surface. 
Thermodynamic calculations indicate the feasibility of these reactions at active pits. Kinetic calculations predict pore depths 
close to those observed. 


L INTRODUCTION 

Reinforced carbon/carbon (RCC) material is a remarkably effective thermal protection material, used in the leading 
edge and nose cone of the space shuttle. This material consists of woven carbon fibers, impregnated with a liquid carbon 
precursor to fill the voids and then converted to carbon by pyrolysis. Such a material exhibits excellent high temperature 
mechanical properties, but as is well known, requires virtually perfect oxidation protection. 

Oxidation protection is attained via several coating steps (ref. 1). First a diffusion coating of silicon-rich silicon 
carbide (SiC) is grown at high temperatures using a pack cementation process. Next a coating of tetraethyl orthosilicate 
(TEOS) is applied via a vacuum impregnation and cured to form silicon dioxide (Si0 2 ) to plug surface cracks and fissures. 
The final coating is a Type A sealant, which is a proprietary slurry coating of sodium silicate, SiC, and other compounds. 
Figure 1 is a micrograph of a polished cross section of this material. The difference in thermal expansion of the carbon 
substrate and SiC coating leads to cracking of the coating and the fluid glass from the TEOS and Type A sealant is 
intended to fill these cracks up to the operating temperature. 

During a shuttle mission, the RCC material is subjected to high temperature aggressive environments (ref. 2). 

Before the launch, the shuttle craft is exposed to the salt laden air of Florida for periods of up to a month. More salt is 
expected to deposit on the wings rather than the nose cone, as the latter is covered by a canopy on the launch pad. These 
salt deposits are likely to remain on the wings during launch and re-entry, especially if they are trapped in surface cracks 
and recesses. Launch lasts only for about 8 min and involves temperatures up to 866 K, temperatures far below the melting 
of boiling points of the salts. Higher temperatures are attained during re-entry, which lasts about 30 min and involves 
temperatures up to 1977 K. The time to reach temperature is about 5 to 8 min and time at the highest temperature lasts 7 to 
8 min. Possible trans-Atlantic (TAL) abort trajectories would involve temperatures to 2095 K. 

The re-entry atmosphere is complex, containing a mixture of atoms and molecules of oxygen, nitrogen, and carbon 
dioxide in a pressure range of about 0.005 to 0.10 atm. These conditions are best simulated with an arc-jet (ref. 1). 
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During a mission, a n umb er of chemical reactions may lead to material recession and/or localized degradation, such 
as the pinholes which have been observed on the wing leading edge material after several missions. The purpose of this 
report is to provide a detailed description of each of the several possible degradation reactions. Specifically, these include: 

1. Oxidation. Passive oxidation may convert the SiC to Si0 2 (s) and CO(g). Under extreme reducing conditions, 
especially in very localized spots, the silica and sodium silicate may no longer be stable. Under such a condition of greatly 
reduced ambient oxygen pressure, solid Si0 2 will not form. Instead, evaporation of the gaseous suboxide, SiO(g), leads to 
rapid material consumption; this is termed "active oxidation". 

2. Reaction of Si0 2 with SiC. This leads to generation of gaseous products and possible bubble formation at the 
SiC/glass interface. 

3. Flow of the sodium silicate. 

4. Vaporization of the sodium silicate layer. 

5. Sodium chloride (NaCl) induced corrosion. This could involve a reaction of both silica and SiC with chlorine 
from the NaCl to form volatile SiCl x (x = 1 to 4) compounds. This process may lead to the observed pinholes, as is 
analyzed in detail in this report. 


H. OXIDATION 

Most metals and silicon-based ceramics (e.g., silicon carbide) are not stable in air at high temperatures and form a 
thin oxide film or scale which acts as a barrier against further attack. There are many reviews of this field (ref. 3 to 5). 
Examples of protective oxides are NiO on Ni, A1 2 0 3 on Ni-Al alloys, and Cr 2 0 3 on Ni-Cr alloys. Some materials do not 
form a solid oxide and hence are unstable in oxygen without a secondary coating. Carbon is, of course, the best example of 
this, since it oxidizes to form CO(g) and C0 2 (g) in a ratio dependent on the temperature and availability of oxygen. 

The degree of protection offered by an oxide film depends on its inherent coverage, adherence, and diffusion 
characteristics. The best protective films have the lowest diffusion coefficient for oxygen; these include alumina (A1 2 0 3 ) 
and silica (Si0 2 ). In high temperature oxidizing gases, SiC forms a thin protective film of silica as: 

SiC(s) (2) + 3/2 0 2 (g) = Si0 2 (s) + CO(g) (1) 

Since oxygen is the mobile species in Si0 2 , both Si0 2 and CO form at the SiC/Si0 2 interface. Figure 2 is a 
schematic illustration of a scale of Si0 2 on SiC, showing the fluxes of 0 2 in and CO out, as occurs during the oxidation 
process. There is general agreement (ref. 5) that the oxidation rate of SiC is limited by the inward diffusion of the 0 2 , so 
the growth rate of Si0 2 is inversely proportional to the film thickness. This results in parabolic scale growth as shown in 
figure 3 for the oxidation of SiC at 1300 °C. Scale growth rates can be measured in either scale thickness or net weight 
gain of oxygen. 

Many factors may influence the diffusivity of oxygen through the silica scale (ref. 5). There is general (but not 
complete!) agreement in the literature that oxygen diffuses through amorphous silica as a molecule and through crystalline 
silica as an oxygen anion, and that the diffusivity of oxygen in amorphous silica is greater than in crystalline silica. 
Amorphous silica tends to form below about 1650 K and tends to crystallize at higher temperatures and/or in the presence 
of impurities. Grain boundaries in crystalline Si0 2 may provide rapid diffusion paths. The form of oxygen also has an 
influence on its diffusion rate and corresponding scale growth rate. Atomic oxygen should diffuse through silica faster than 
molecular oxygen, leading to more rapid scale growth rates. 

An important consequence of a diffusion-controlled growth of a thin product scale is a steep gradient in oxygen 
potential through the scale. For silicon carbide oxidizing in one atmosphere of oxygen, the oxygen pressure at the gas/oxide 
interface is one atm, and the oxygen partial pressure at the oxide/carbide interface can be as low as 4.3x10” atm at 
1500 K. This is a significant issue in determining the reactivity at this interface. 

The protection scheme for the RCC material is unique in that a film of Si0 2 is applied via a vacuum impregnation 
process to the SiC. At high temperatures, further oxidation of the SiC particles in the Type A sealant and silicon-rich SiC 
coating on the carbon/carbon is likely to occur via reaction (1). Parabolic oxidation kinetics require that additional silica 
formation slows the rate as the scale thickens. The oxidation reaction also generates CO(g) as a product, which either 
dissolves and diffuses through the glass or else forms bubbles. Over very long times and at high temperatures, CO 
evolution may cause scale rupture and lead to oxidative consumption of the SiC coating. 


2 



At very low oxygen pressures, e.g., beneath a stagnant boundary layer of an inert gas, the passive Si0 2 (s) becomes 
unstable and active oxidation may consume the SiC. The transition pressure is derived in Appendix A. Note that active 
oxidation only occurs for extreme conditions of high temperature and reduced oxygen availability. 

In summary, continued thermal oxidation adds silica to the glass coating on the RCC and therefore slows the 
oxidation rate. Under very low oxygen partial pressures and high temperatures, however, active oxidation and rapid 
material consumption can occur. 


HI. REACTION OF SiC AND Si0 2 — POSSIBLE BUBBLE FORMATION BENEATH THE SCALE 

The silicon carbide/silica system is somewhat unique in that at high temperatures, these components react according 
to equation (2) to generate volatile products at their mutual interface: 

SiC + 2 Si0 2 = 3 SiO(g) + CO(g) (2) 

At very high temperatures, this reaction can generate high gas pressures at the SiC/Si0 2 interface (refs. 7 and 8). If 
the gases do not escape by diffusion, they may form bubbles and disrupt the Si0 2 -glass coating. Generally it is assumed 
that when the total pressure of SiO(g) plus CO(g) exceeds somewhat the ambient pressure, a gas bubble will form. 

However, the exact conditions for gas bubble formation are much more complex, involving surface tension, nucleation 
sites, viscous flow, and other factors. The pressures and relative amounts of SiO(g) and CO(g) generated by the above 
reaction depend on temperature and the thermodynamic activities of Si and C in SiC; in general the 3:1 ratio for SiO to CO 
gases of equation (2) is not satisfied. 

The exact value for the carbon activity in SiC can be calculated from the standard Gibbs energy of formation, 
AG° SiC and the corresponding equilibrium constant, K: 

Si + C = SiC AG° = -61.585 kJ/mol at 1500 K = -RT In K 3 (3) 

K 3 = 139.53 = l/fasjac) (4) 

For a Si = 1 ac = 7.167xl0 -3 

Thus the carbon activity can only vary between 1 (carbon-saturated SiC) and 0.007167 (silicon-saturated SiC) at 
1500 K. Correspondingly, the silicon activity can only vary between 0.007167 (carbon-saturated SiC) and 1 (silicon- 
saturated SiC). If an activity of one component is set, the conresponding activity of the other component is decided. For the 
reaction of SiC and Si0 2 , the pressures of SiO(g) and CO(g) can be calculated for both silicon- and carbon-saturated SiC. 
The best way to illustrate this is to examine the SiO and CO isobars over the Si-C-0 isothermal stability diagram. A recent 
article by Gaskell (ref. 9) describes the construction of these diagrams in detail. 

Stability diagrams with log (ac) and log (P Q ) as axes are the easiest to construct and probably the most relevant 
for the RCC case. The diagram shown in figure 4 presents a map of the stable condensed phases for a given carbon activity 
and oxygen pressure. We are interested in the SiC/Si0 2 interface, which is represented by the line AB in the stability 
dia gram Any point on this line can be specified by setting an activity of carbon or a partial pressure of oxygen. Point A is 
carbon-saturated SiC and point B is silicon-saturated SiC. 

In practice, identifying a carbon activity and a partial pressure of oxygen by direct measurement is difficult. So we 
construct isobars of SiO(g) and CO(g) on figure 4. From the isobars, we can approximate the pressures of SiO and CO for 
carbon-saturated SiC, silicon-saturated SiC, as well as intermediate compositions of SiC. These pressures are given by the 
intersections of the SiO(g) and CO(g) isobars at the appropriate point on line AB. These values can be calculated for points 
A and B, respectively, from the following expressions (ref. 7): 

Carbon-saturated SiC: 

For P sio : SiC + Si0 2 = C + 2 SiO(g) K 5 = P Si0 2 
For P c0 : 3C + Si0 2 = SiC + 2 CO = P qq 2 
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Silicon-saturated SiC: 


For P Si0 : Si + Si0 2 = 2 SiO(g) K ? = P Si0 2 (7) 

For P co : 2 SiC + Si0 2 = 3 Si + 2 CO K g = P co 2 (8) 

Table I lists these pressures for a range of temperatures. These results are shown graphically in figure 5, which is a 
plot of total pressure of CO(g) and SiO(g) for silicon-saturated SiC + Si0 2 , carbon-saturated SiC + Si0 2 , and SiC + Si0 2 , 
which produces SiO(g) and CO(g) in a 3 to 1 molar ratio according to reaction (2). For figure 5, one atmosphere has been 
taken as a point of bubble formation. Note that P(SiO) + P(CO) is significantly higher for carbon-saturated SiC than for 
silicon-saturated SiC. This important difference amounts to about a 250 K temperature advantage for silicon-saturated SiC 
over carbon-saturated SiC. 

Silicon-saturated SiC reacts much less with the Si0 2 glass to generate gaseous products than does carbon- 
saturated SiC. For this reason the SiC diffusion coating on the carbon/carbon is silicon-rich. This is also the reason 
why Si particles would be preferred to SiC particles in any slurry sealant intended to plug holes in the SiC coating. 


IV. FLOW OF THE GLASS 

The Type A sealant is a sodium silicate glass. There is clearly a correspondence of viscosity to temperature and 
certainly to Na content as well. At temperatures below the maximum re-entry temperature, the glass is molten and tends to 
flow. This is beneficial in that it will fill any cracks in the SiC that result from the thermal expansion mismatch with the 
RCC. However, if NaCl is present at the base of these cracks, the flow of glass will not displace this NaCl which may still 
lead to some damage. This will be discussed in section VI. In addition, if the glass flows too much, it may be blown off 
the surface by viscous drag, exposing the SiC and possibly the carbon/carbon itself to attack. 


V. VAPORIZATION OF THE GLASS 

Let us first consider Si0 2 and the gas phase alone. Figure 6 shows a volatility diagram for Si0 2 on Si at 1800 K. 
We are concerned with the Si0 2 portion of this diagram as it exists on an SiC substrate or coating. Several references 
describe the construction of these diagrams (ref. 4, 10 and 11). For a given temperature, such a diagram plots the partial 
pressures of the various vapor species above a condensed phase as a function of oxygen pressure. Note that Si0 2 (s) does 
not only vaporize to Si0 2 (v), but rather to several different species. Increasing the oxygen pressure suppresses vaporization 
of SiO, but has no effect on the vaporization of Si0 2 to Si0 2 (g). Qualitatively, this can be seen from the reactions: 


Si0 2 (s) = SiO(g) + 1/2 0 2 (g) 

(9) 

with d log P Si0 /d log Pq 2 = -1/2 


Si0 2 (s) = Si0 2 (g) 

(10) 

with d log P Si o 2 /d log Pq 2 = 0 



The outer coating on the RCC materials is the Type A sealant, which is primarily a sodium silicate glass. There are 
several sodium silicate glass compositions with different ratios of Na^ to Si0 2 . Since this glass is in contact with the Si0 2 
from TEOS, it likely becomes the most silica-rich composition — Na^^Oj. We can construct a volatility diagram for 
NajSijOj. However, this is a more complex system than Si0 2 and can be handled another way. 

In the case of a sodium silicate glass, Na is expected to vaporize preferentially. Since this involves a decomposition, 
oxygen is expected to suppress the volatility of Na The best way to calculate vapor pressures is with a free energy 
minimiz ation type computer code (refs. 12 and 13). Such a code considers all the atoms from a given system and 
determines the products with the minimum free energy. The vapor pressure can be calculated for a vacuum or for a system 
with a given over-pressure of gas. Figure 7 shows that calculated vapor pressures above Na^jS^Oj in vacuum and figure 8 
shows the calculated vapor pressures above Na 2 Si 2 O s in 0.05 atm oxygen. Note that Na is the primary vapor species and 
that 0.05 atm oxygen suppresses the vapor pressure of Na by more than two orders of magnitude. These are thermodynamic 
(equilibrium) vapor pressures. 


4 



From the standpoint of kinetics, the maximum possible vapor flux (ref. 14) for a given vapor pressure can be 
calculated from the Hertz-Langmuir expression: 


J = aP/(2nMRT) 1/2 (11) 

where J is a flux in mole/(unit area-unit time), a accounts for stearic and other kinetic limitations and any backflux of 
species (ref. 14), M is the molecular weight of the gas, R is the gas constant and T is the absolute temperature. The total 
evaporation rate comprises the sum of fluxes calculated for each vapor species. This maximum possible vapor flux of 
equation (11) is observed in a vacuum or rapidly flowing gas with no boundary layer. In most problems, a boundary layer 
is present and the evaporation is limited by gaseous diffusion, as described by the expression: 


J = hJP-Pg) (12) 

Here P is the pressure of the vaporizing species in the carrier or ambient gas (which is often negligible) and is 
the mass transfer coefficient. Semi-empirical values for h,,, are tabulated for a variety of geometries (ref. 15). 

The experimental measurement of vapor pressures and evaporation fluxes is a large area of research (ref. 16). The 
easiest way to measure fluxes is with interrupted or continuous weight loss. Continuous weight loss involves suspending the 
specimen from a balance in a furnace and continuously monitoring the weight change; a schematic of such a system is 
shown in figure 9. In these types of measurements, there are two approaches. One involves measuring "Langmuir" 
evaporation from a free surface into a vacuum. The other involves containing the sample in a Knudsen cell, which is a 
small enclosure with a well defined orifice for vapor sampling. This enclosure allows the gas and the condensed phase to 
approach equilibrium, so that a near equilibrium vapor pressure evolves from the cell. 

Because weight-change measurements only give total vapor fluxes, some method of analysis is necessary to 
determine the vapor composition. Perhaps the most valuable method is mass spectrometry (ref. 14). A Knudsen cell is 
generally used in these systems, since the vapor must be columnated as a molecular beam. Such a mass spectrometer can 
measure pressures only as high as about 10~^ atm. Sampling higher pressures requires a series of differentially pumped 
vacuum chambers (ref. 15). These mass spectrometric systems can sample condensible gases, such as SiO, Na, and NaO as 
well as noncondensible gases, such as CO, C0 2 , and 0 2 . A gas chromatograph can also be used at higher pressures, but it 
can only sample noncondensible gases. 

Another method of studying vaporization involves rapid laser heating and using a mass spectrometer to probe the 
imposition of the resultant plume (ref. 18). Such a test may be more relevant to the rapid heating encountered in re-entry. 

In summary, the Type A sealant is expected to vaporize preferentially to lose Na. The pressure of oxygen in 
the re-entry environment will suppress this to some degree. Vapor pressures can be calculated with free energy 
minimization computer codes and the corresponding maximum possible weight loss rates are calculated from the 
highest vapor pressures. Vapor fluxes and vapor compositions can be measured experimentally using a gravimetric 
balance and a mass spectrometer, respectively. 


VL NaCI INDUCED CORROSION 

It is well-known that NaCI is an aggresive species in a variety of corrosive environments. In low temperature 
aqueous environments Cl' from NaCI is an aggressive anion, promoting pitting (refs. 19 and 20). In high temperature 
combustion environments, NaCI combines with sulfur impurities in the fuel and water vapor to form the highly stable 
Na2S0 4 molecule. Deposits of fused Na2S0 4 can be highly corrosive (refs. 4 and 5). 

In the case of the RCC material on the shuttle, evidence suggests that yet another form of NaCl-induced attack may 
be occurring. As mentioned, the shuttle can reside on the lauch pad for weeks with the RCC material of the wings exposed 
to the moist salty seabreezes. Thermal excursions of the RCC material lead to extensive craze cracking of the solidified 
Type A sealant. A micrograph of this is shown in figure 10. Generally, these surface flaws are not repaired or filled 
between flights, and therefore they can serve as sites where sea salt is retained and protected from loss by the aerodynamic 
forces upon ascent. In the analysis which follows, sea salt stored in surface cracks is considered to be heated upon reentry 
to high temperature where degradative reactions can result. Pinhole formation is suggested to result from a cyclic 
chloridation/oxidation of silicon carbide, with the formation of a condensed nonprotective silica fume. 
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VI A. THERMODYNAMICS 


As discussed in section V, the computer-assisted analysis of multicomponent, multiphase equilibria between several 
condensed phases and a gas phase can be performed by a free energy minimization program to calculate the partial 
pressures of the product vapor species (refs. 12 and 13). To some extent, the results depend upon the phases that are 
assumed to be present, but the amounts of stable condensed phases do not matter, so long as the phases are retained at unit 
thermodynamic activity. Likewise, the partial pressures of the nitrogen and oxygen gases and the total gas pressure should 
not be important to this thermodynamic analysis, because volatile nitrogen compounds are not very stable, and the oxygen 
is already accounted for by assuming that silica is present (or not present). For the calculated results presented in tables II 
and VII, covering the temperature range from 1200 to 1700 °C (2192 to 3092 °F), the condensed phases SiC and NaCl are 
gggnmpH to be present, in equilibrium with either carbon or else with silicon, and either in the presence of silica (passive 
reaction) or else in the absence of silica (active reaction). The ambient gas atmosphere is assumed to be pure nitrogen at 
1 atm pressure, because some specification of total pressure is required by the program. In fact, for the thermodynamic 
considerations, the total pressure and the ambient gas composition do not affect the calculated results if the condensed 
phases are specified. 

From the Si-C-0 phase stability diagram of figure 11 for 2000 K, the coexistence of SiC and Si0 2 can only lie 
along a line limited by a further equilibrium with either carbon or else with silicon. Then to examine the full range of 
values possible, tables II to VII provide the results of calculations for the vapor pressures of species for the 4-phase 
equilibrium of NaCl + SiC + Si0 2 + C (point A of fig. 12), or else for NaCl + SiC + Si0 2 + Si (point B). In tins case, the 
microenvironment in the pinhole is assumed to be sufficiently oxidizing that Si0 2 exists (passive reaction). Additionally, 
for the proposed case that the local pinhole environment actually becomes "active" due to the total evaporative loss of 
Si0 2 , the 3-phase equilibria for NaCl + SiC + C (point D) and for NaCl + SiC + Si (point E) have also been evaluated. 
These four sets of equilibria should bracket any possible set of service conditions for passive oxidation, and for active 
oxidation, in the NaCl-contaminated pinhole. 

Figure 12 su mmar izes the calculated results for the points A and B of figure 11, as presented in detail in tables II 
and VII for the condition of passive reaction, i.e., that silica in the pinhole is not totally consumed by the evaporation 
reactions. For any given temperature, the highest vapor pressures are exhibited by NaCl, which species can boil out of the 
crack (depending upon the ambient temperature and pressure). The partial pressures for CO, SiO, SiCl 2 , etc. are plotted as 
a small range of values for a given temperature, corresponding to the total possible range between equilibrium with carbon 
or else with silicon. For a qualitative understanding of the degradation mechanism, this range for the values is not 
important. However, as the temperature is raised, the magnitudes of both SiO and SiCl 2 vapor are high enough to support 
a significant flux of silicon from the base of the pit to the surface. 

Figure 13 su mmarizes the detailed results for the points D and E of figure 11, as presented in detail in tables II to 
VH, for the condition of active (silica-free) reaction in the pinhole. Again the range in the possible values is not very 
important. In form and magnitude, the results for figure 13 are rather similar to those for figure 12. Thus, in the presence 
of condensed NaCl, except for the absence of the volatile oxides SiO and CO, a large difference is not expected between 
"passive" and "active" vaporization of SiC in a pinhole. 


VI B. CIILORIDATION/OXIDATION REACTION MECHANISM 

Figure 14 shows schematically a reaction sequence which rationalizes the generation and growth of pinholes 
through the SiC coating on carbon/carbon. In figure 14(a), sea salt is condensed into the cracks and crazes of the SiC 
coating as the shuttle is exposed to the ambient Florida environment. Since the cracking and salt deposition occurs at low 
temperature, only very little oxide should be present in the crack. During the minor heating upon ascent, in the face of high 
aerodynamic pressure, the salt would be retained in the crack until reentry, when the high temperature would melt the salt 
and initiate the reactions indicated in figure 14(b), leading to an evaporative loss of NaCl, SiO, CO and SiCl 2 . This 
transient reaction stage corresponds to the "passive" reaction of SiC/Si0 2 described in figure 14(b). Because of the very 
limited amount of oxygen in the crack, all this oxygen should be lost by SiO and CO formation before the NaCl is lost. 
Thereafter, as shown schematically in figure 14(c), the base of the pit should suffer ' active chloridation leading to SiCl 2 (v) 
loss. 

At the external surface, the SiCl 2 vapor would encounter oxygen by collisions with the ambient gas phase. Because 
in contact with some oxygen, condensed silica is a much more stable phase, the volatile SiCl 2 (and any SiO) species would 
react to form a solid (particulate) fume of Si0 2 , as shown schematically in figure 14(c). The relevant reaction at the 
external surface is: 
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(13) 


SiCl 2 (v) + 0 2 (g) = Si0 2 (fume) + 2 Cl(g) 

The atomic chlorine product species are trapped within the pinhole, and diffuse back to the base of the pinhole to 
repeat the reaction to form SiCl 2 (v), as shown schematically also in figure 15(a). There are several possible mechanisms for 
the removal of carbon from SiC at the base of the pinhole. The carbon could diffuse away into the Si-saturated SiC, or else 
it might be removed as a volatile species, e.g. some hydrocarbon molecule if any hydrogen is present in the pit In any 
case, carbon removal is not thought to be involved in rate control. The fume particles at the surface constitute a porous, 
non-protective deposit. (Actually, such condensed fume is seen smeared across the external surface of the SiC coating of 
the shuttle leading edge) The formation of this porous surface deposit is important to the mechanism, because it limits the 
entry of oxygen from the ambient into the pinhole, an action which would otherwise repassivate the SiC by silica formation 

in the pinhole. 

Such a chloridation/oxidation reaction model has many analogs in systems that have been well studied. The 
generation of fume resulting from the evaporation of liquid metals was described expertly by Turkdogan, Grieveson and 
Darken (ref. 21). Hinze and Graham (ref. 22) found copious Si0 2 fume as a product of the active oxidation of solid Si. 
Internal vapor transport cycles are well known in the chloridation of alloys (refs. 23 to 26). If one excludes the electro- 
chemical aspects for materials attack in electrolytic media, this mechanism is quite similar to the localized aqueous 
corrosion (pitting) of alloys which form passive films in oxygenated environments (refs. 19 and 20). During the pitting of 
such alloys, chloride and hydrogen ions are concentrated in the pits, where the local environment differs greatly from the 
bulk solution, and the alloy suffers active corrosion. 


VI C. KINETICS OF THE CHLORIDATlON/OXIDAnON REACTION-RATE OF PINHOLE GROWTH 

To analyze the rate at which pinholes should grow/deepen according to the mechanism just proposed, one needs to 
decide whether the diffusion of reactant and product gas species in the pinhole occurs by Knudsen (free molecular) 
diffusion (for L » d) or else by viscous diffusion (for \ « d), where \ is the mean free path for species i and d is the 
characteristic dimension of the system (diameter of the pinhole). The following parameters bracket the conditions in the 
pinhole: 

Temperature range: 1500 to 2100 K 

Diameter range of pinholes: 380 to 510 mm (0.015 to 0.020 in.) 

Total pressure in pinhole: 0.01 to 0.07 atm 

Actually, the pressure in the pinhole has been estimated on the basis of the shock wave against the wing leading edge. 


VI C.l. MEAN FREE PATHS 

To calculate the mean free path of a species in the gas phase, the diameter of the species is required, but these are 
not always available. From reference 27, the diameter of N 2 is 3.798 Angstroms. Although there are no data available for 
the diameter of SiCl 2 (v), its diameter can be estimated by a linear interpolation between SiCl(v) and SiCl 4 (v). Reference 28 
shows that the diameters for CC1 X , BC1 X , and CF X are linearly dependent upon x. According to this interpolation, the 
molecular diameter for SiCl 2 (v) is 4.491 Angstroms, as listed in table VIII. 

The mean free path of a species i at temperature T is given by (ref. 28): 

1; = l/[(2) 1/2 7tn i d i 2 ] = k B T/[(2) 1/2 JtP i d i 2 ] (14) 

where n- and P are the density and partial pressure of i, and k B the Boltzmann constant. For the model and calculations to 
follow, the interdiffusion of the vapor species N 2 and SiCl 2 are considered to be most important. To consider collisions 
between N 2 and SiCl 2 , the summation of the radii for N 2 and SiCl 2 is set as d ; and P N ,, as Pj in equation (14); the 

calculated results are listed in table IX. . 

The mean free paths of SiCl 2 should be somewhat smaller than those listed in table IX because of the collisions 

between SiCL molecules. According to the plot of partial pressure profiles of figure 15(b), atomic chlorine has a range of 
pressure which is nearly twice that for SiCl 2 , but its atomic diameter is smaller (table VIII), and its diffiisional transport 
should not be rate limiting. The measured diameters of the pinholes are in a range of 380 to 510 mm, which are much 
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larger than the mean free paths for SiCl 2 at temperatures of 1500 to 2100 K. Therefore, viscous diffusion of SiCl 2 should 
control the transport of SiCl 2 in pinholes and not Knudsen (free molecular) diffusion. A simplification of the microenviron- 
ment of the pinhole gas phase to a N 2 -SiCl 2 binary system can be made for calculation of the diffusivity of SiCl 2 in a 
pinhole. 


VI C.2. DIFFUSION CALCULATIONS 

The model shown in figure 15(a) was proposed for the growth of a pinhole in SiC. The gradient for SiCl 2 pressure 
in the pinhole supports the diffusion of SiCl 2 through the N 2 "background gas" outward to the surface, as shown in 
figure 15(b). Equation (13) describes the fume-forming reaction at the external surface. Theories for calculating the 
diffusivity in a binary system are still not well developed, but some empirical equations are available for calculation of the 
diffusivity. The following equation is used to calculate the diffusivity of SiCl 2 in N 2 (ref. 15): 


( M/2 


D SiCl 2 -N 2 = — | 


v* 3 y 


M/2 


2 m v; , 


■SiCU 


2m 


% 


- 3/2 


r d SiCl 2 + d N^ 


(15) 


Here D SiC1 _ N is the diffusivity in cm 2 /sec, k B is the Boltzmann constant, ntj is the weight per molecule (i.e., 
molecular weight/Avogadro Number), T is the absolute temperature, p is the total pressure, and d ; are the molecular 
diameters. The calculated results are listed in table X. The flux of SiCl 2 (v) in the pinhole shown in figure 15 is given by: 


dC SiCl 2 

J SiCl 2 - D SiCl 2 —%— 


D SiCl 2 P SiCl 2 

pFF 


(16) 


where D SiC1 is the diffusivity calculated in equation (15), dC sic , 2 /d£ is the concentration gradient, £ is the depth of the 
pinhole, anc^P SiC , 2 ; s ^ vapor pressure of SiCl 2 in the presence of NaCl and SiC, but in the absence of 0 2 . The pressure 
of SiCl 2 is from carbon-saturated SiC, since the chloride forming reaction releases carbon (fig. 14(c)). Every mole of SiCl 2 
formed corresponds to one mole of SiC removed, so J SiCl2 = J SiC . The growth (penetration) rate of the pinhole is related 
to the flux as follows: 

^ - v / v DsiCl 2 PsiCl 2 (17) 

Tt ' V SiC J SiC' v SiC 


where V siC is the molar volume of SiC (12.43 cm 3 /mol). Upon integration of equation (17): 


2 2 V SiC D SiCl 2 P SiCl 2 




RT 


t • k p t 


(18) 


The parabolic rate constants k^ for the pinhole growth are calculated and listed in table XI. The calculated pinhole growth 
kinetics are plotted in figure 16 for 1200 to 1700 °C at total pressures of 0.01 and 0.07 atm. The extreme condition for the 
shuttle service can amount to 7.5 to 8 min at a temperature of 1600 °C. These conditions have been highlighted in 
figure 16, and they correspond to a calculated depth for pinhole penetration of 490 to 1340 pm (19 to 53 mils). Three 
specific pinholes have been studied in some detail (ref. 29), two pinholes in R/H panel #7 had depths of 15 and 28 mils. 
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and a pinhole in R/H T-seal R/S 11 had a depth of 15 mils. Considering the uncertainties and approximations, such 
agreement between analysis and observation may be considered fortunate. On the other hand introduction of a tortuousity 
factor into equations (17) and (18) to account for the turns in the capillary brings about very good agreement. Indeed, the 
chloridation/oxidation mechanism represents a logical interpretation for pinhole formation. Clearly, the presence of a fume 
at the pinhole entrance/exit has been observed at the RCC surface of the shuttle. 

In summary, it is possible that the observed pinholes form from the following mechanism. The sea salt 
deposited on the exposed wing surfaces very likely collects in the cracks in the glass sealant This can penetrate to 
the SiC surface and react to form silicon chlorides during re-entry. The silicon chlorides are oxidized to Si0 2 fume 
at the surface and the chlorine is released to form more silicon chlorides. Thermodynamic calculations indicate the 
feasiblity of such reactions. Kinetic calculations predict pit depths close to those measured. 


VI C.3 STOPPAGE AND AVOIDANCE OF PINHOLES 

The chloridation/oxidation mechanism presented here provides answers for a number of problem questions which 
are tied to pinholes in the SiC coating on the shuttle: 

Question 1: Why have p inho les never been observed or simulated in the arc jet or other thermal exposure testing 
of the materials? Answer: Such test coupons were not provided with imbedded seasalt prior to exposure. 

Question 2: Why does the nose cone, which experiences the same thermal service cycle as the leading wing 
edges, suffer only limited pinhole formation? Answer: On the launch pad, the nose cone is covered by a canopy which 
prevents contact with salt spray. 

Question 3: How can one prevent pinhole formation? Answer: As soon as possible after the shuttle lands, the 
wing panels should be sprayed or painted with some sealant to fill in the surface cracks so toat salt contamination is not 
gathered into surface cracks either in transport or in position on the launch pad. A vacuum impregnation is recommended. 
A canopy should cover the wings of the shuttle on the launch pad until shortly before takeoff. 

Question 4: What can be done to clean out existing pinholes and prevent further growth? Answer: Perhaps one 
can clean out the existing pinholes with an element which forms a quite stable volatile chloride, oxide, and oxychloride so 
that any residual chlorine can be extracted from the pinholes without the formation of a blocking porous oxide fume 
deposit. The elements Cr and Mo both satisfy this requirment, although Cr is preferred. If the open pinholes were provided 
with a chromic acidic wash, then upon heating in the presence of oxygen, Cr0 3 (v) and Cr0 2 Cl 2 (v) should form to clean 
out the pinhole at temperatures lower than those where SiCl 2 (v) reaches a high vapor pressure. 


VH. SUMMARY AND CONCLUSIONS 

The RCC mater ial is exposed to high temperatures and aggressive environments at several points in a typical 
mission. Several possible degradation modes have been discussed: 

1. Passive oxidation of the SiC coating and SiC particles in the slurry. At extreme conditions of high temperature 
and low oxygen pressure, the transition to rapid active oxidation may occur. 

2. Reaction between SiC and Si0 2 . For carbon-saturated SiC this reaction would lead to a high CO(g) pressure 
generated at the SiC/Si0 2 interface; silicon-saturated SiC is best for minimizing this pressure. 

3. Vaporization of the glass. The sodium component is expected to vaporize preferentially. This would be 

minimized by an overpressure of oxygen. . 

4. NaCl induced corrosion. This may act to volatilize the protective glass and lead to local pitting. A 
chloridation/oxidation mechanism causing active reaction of SiC at the base of the pit leads to the deposition of a porous 
Si0 2 fume at the mouth of the pit and recycling of the atomic chlorine reactant. Extended or repeated operation of this 
degradation mechanism can explain the pinholes observed in the shuttle leading wing surfaces. 
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TABLE I.— PARTIAL PRESSURE (atm) RANGE FOR SiO AND CO AT 
THE SiC/VTTREOUS Si0 2 INTERFACE 


Temperature 

K(°F) 

p Sicr carbon 

saturated 

Pqj - carbon 
saturated 

Psio-stlicon 

saturated 

P co -carbon 

saturated 

1200 (1701) 

2.23xl0~ 8 

6.07xl0" 5 

5.47xl0" 7 

4-lOxlO" 9 

1300 (1881) 

4.13xl0~ 7 

6.08x10^ 


9.54x1 0" 8 

1400 (2061) 

5-OOxlO" 6 


7.28x1 0 -5 


1500 (2241) 

4.30x1 0 -5 

2.37X10 -2 

5.09X10 -4 

1.43xl0" 5 

1600 (2421) 

2.81 xlO -4 

1.04x10"’ 


1.09X10" 4 

1700 (2601) 


3.82x10"’ 

|gggjj| 

6.79X10" 4 

1800 (2781) 

6.24xl0" 3 

1.21 


4.42xl0" 3 

1900 (2961) 

2.28xl0" 2 

3.35 

1.19x10"’ 

2.35xl0" 2 

2000 (3141) 

7.29xl0~ 2 

8.34 

3.14x10"’ 
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TABLE 0.— CALCULATED ACTIVITIES FOR GASEOUS AND 
CONDENSED SPECIES WITH DIFFERENT CONDENSED 
PHASES PRESENT AT 1200 °C (2192 °F) 


Species 



SiC + NaCl + Si 

SiC + NaCl + C 

N 2 

0.906752 

0.89744 

0.90704 

0.907346 

NaCl 

0.057474 

0.05747 

0.05747 

0.0574737 


0.035083 

0.03508 

0.03508 

0.0350828 

SiO 

0.000275 

0.00988 

0 

0 

Na 

0.000269 

6.60x1 0" 5 

0.00027 

6.61 xl0~ 5 

SiCI, 

0.00013 

3.20xl0‘ 5 

0.00013 

2.84xl0 -5 

CO 

1.49xl0~ 5 

2.80x1 0' 5 

0 

0 

SiCl 3 

2.50x10"* 

2.20x10-* 

2.50x10-* 

2.22X10 -6 

SiCL 

1.65xl0' 7 

5.90xl0 -7 

1.60X10 -7 

5.94xl0 -7 

Cl 

6.54xl0" 8 

2.70xl0 -7 

6.50xl0" 8 

2.66xl0 -7 

SiCl 

3.24xI0 -8 

8.40xl0 -8 

3.20xl0 -8 

1.73xl0' 9 

Si 

5.17xl0 -9 

1.70X10 -9 

5.20xl0“ 9 

6.80x10-" 

Naj 

2.85xl0 -9 

1.70xl0 -1 ° 

2.90X10 -9 

1.72X10 -10 

co 2 

1.45xl0" u 

6.80x10'" 

0 

0 

SiO, 

3.21 xlO” 12 

3.10x10'" 

0 

0 

Cl 2 

1.84xl0 -12 

3.20xl0 -12 

1.80xl0 -12 

3.05x10-" 

NO 

8.53xl0' 15 

7.40xl0 -14 

0 

0 

NaO 

2.08xl0 -16 

4.50xl0' 16 

0 

0 

SiC 

4.60xl0~ 18 

4.60x1 0 -1 8 

4.60xl0 -18 

4.60xl0” 18 

n 2 o 

4.97xl0 -19 

4.29xl0' 18 

0 

0 

0 2 

9.79xI0 -24 

7.44xI0~ 22 

0 

0 

no 2 

3.21xl0“ 28 

2.43x1 O' 2 * 

0 

0 

Condensed 





SiC 

1 

1 

1 

1 

Si0 2 

1 

1 

0 

0 

Si 

1 

0.01315 

1 

0.0131503 

NaCl 

1 

1 

1 

1 

Na 2 Si 2 O s 

0.000197 

0.0001 

0 

0 

Na 2 Si0 3 

2.4x1 0' 5 

1.20xl0 -5 

0 

0 

C 

0.01315 

1 

0.01315 

1 
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TABLE m.— CALCULATED ACTIVITIES FOR GASEOUS AND 
CONDENSED SPECIES WITH DIFFERENT CONDENSED 
PHASES PRESENT AT 1300 °C (2372 °F) 


Species 


SiC + NaCl + 
Si0 2 + C 

SiC + NaCl + Si 

SiC + NaCl + C 

n 2 

0.777301 

0.733075 

0.779007 

0.780047 

NaCl 

0.144609 

0.144609 

0.144609 

0.144609 

Na^CI-j 

0.0749444 

0.0749444 

0.0749444 

0.0749444 

SiO 

0.00157338 

0.00022283 

0 

0 

Na 

0.00096278 

0.00026971 

0.00096278 

0.00026971 

SiCl, 

0.00046709 

0.00011938 

0.00046709 

0.00011938 

CO 

0.0001328 

0.0467481 

0 

0 

SiCl 3 

8.76x10-* 

7.99x10“* 

8.76x10-* 

7.99x10 -* 

SiCL 

4.31xl0' 7 

1.41X10 -6 

4.31X10” 7 

1.41X10 -6 

Cl 

3.75xl0 -7 

1.34x10-* 

3.75xl0” 7 

1.34x10-* 

SiCl 

2.60X10' 7 

1.86x1 O' 8 

2.60x1 0" 7 

1.86xl0" 8 

Si 

5.23xl0~ 8 

1.05xl0' 9 

5.23xl0' 8 

1.05xl0' 9 

Na2 

2.40xl0~ 8 

1.88xl0' 9 

2.40x10-* 

1.88xl0' 9 

co 2 

3.21xl0~’° 

7.99x1 O' 7 

0 

0 


6.33x10'” 

6.33x10'” 

0 

0 

Cl, 

1.64x10'” 

2.09xl0' 10 

1.64x10'” 

2.09xl0' 10 

NO 

1.34xl0" 13 

9.22xl0' 13 

0 

0 

NaO 

7.72xl0' 15 

1.53xl0' 14 

0 

0 

SiC 

2.60xl0' 16 

2.60xl0' 16 

2.60X10' 1 * 

2.60x10'’* 

N,0 

7.09xl0' 18 

4.72x1 O' 17 

0 

0 

°2 

l.llxlO' 21 

5.52xlO' 20 

0 

0 

no 2 

3.99xl0' 2 * 

1.93X10' 24 

0 

0 

Condensed 





SiC 

1 

1 

1 

1 

Si0 2 

1 

1 

0 

0 

Si 

1 

0.0200578 

1 

0.0200578 

NaCl 

1 

1 

1 

1 

Na 2 Si ? 0 5 

0.00050626 

0.00028053 

0 

0 

Na^SiO, 

5.86xl0' 5 

3.25xlO' 5 

0 

0 

C 

0.0200578 

1 

0.0200578 

1 
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TABLE IV.— CALCULATED ACTIVITIES FOR GASEOUS AND 
CONDENSED SPECIES WITH DIFFERENT CONDENSED 
PHASES PRESENT AT 1400 °C (2552 °F) 


Species 

1 


SiC + NaCl + Si 

SiC + NaCl + C 

N 2 

0.558418 

0.456632 

0.563253 

0.564219 

NaCl 

0.302584 

0.2637 

0.304386 

0.305491 

Na,Cl, 

0.126453 

0.0960413 

0.127964 

0.1288954 

SiO 

0.00722118 

0.00125192 

0 

0 

Na 

0.00292722 

0.00084934 

0.00293889 

0.00093826 

SiCl, 

0.00142103 

0.00038532 

0.00142662 

0.00042375 

CO 

0.0009437 

0.181105 

0 

0 

SiCl 3 

2.60x10"® 

2.12xl0 -5 

2.62x1 0" 5 

2.45x10"® 

Cl 

1.74x10"® 

5.23x10"® 

1.75x10"® 

5.49x10"® 

SiCl 

1.62x10"® 

1.46xl0" 7 

1.62x10"® 

1.53xl0" 7 

SiCL 

9.94xl0~ 7 

2.43x10"® 

1.00x10*® 

2.94x10"® 

Si 

4-OOxlO" 7 

1.20X10" 8 

4.00xl0" 7 

1.20xl0" 8 

Na2 

l.53xlO" 7 

1.29xl0" 8 

1.54xl0" 7 

1.57x10"® 

co 2 

5.12X10" 9 

5.67x10"® 

0 

0 

Si0 2 

8.71xl0' 10 

8.71 xlO" 10 

0 

0 

Cl 2 

1.12xlO" 10 

1.01 xlO" 9 

1.12xlO" 10 

l.llxlO" 9 

NO 

1.39xl0" 12 

7.24xl0" 12 

0 

0 

NaO 

1.79xl0" 13 

3.00xl0 -13 

0 

0 

SiC 

9.04x10"'® 

9.04X10" 15 

9.04xl0" 15 

9.04xl0" 15 

N 2 0 

6-lOxlO -17 

2.88xl0" 16 

0 

0 

°2 

7.20xl0" 20 

2.40xl0" 18 

0 

0 

no 2 

2.56x1 0 -24 

7.69xl0" 23 

0 

0 

Condensed 





SiC 

1 

1 

1 

I 

Si 

1 

0.0300563 

1 

0.0300563 

Si02 

1 

1 

0 

0 

NaCl 

0.948752 

0.826832 

0.954403 

0.957867 

Na 2 Si 2 0 5 

0.00109306 

0.0005308 

0 

0 

Na,SiO, 

0.00013177 

6.40x1 0" s 

0 

0 

C 

0.03000563 

1 

0.0300563 

1 
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TABLE V.— CALCULATED ACTIVITIES FOR GASEOUS AND 
CONDENSED SPECIES WITH DIFFERENT CONDENSED 
PHASES PRESENT AT 1500 °C (2732 °F) 


Species 


SiC + NaCl + 
Si0 2 + C 

SiC + NaCl + Si 

SiC + NaCl + C 

N 2 

0.520544 

0.389022 

0.539248 

0.540842 

NaCl 

0.360333 

0.288156 

0.370286 

0.373165 

NazClj 

0.0771595 

0.0493442 

0.0814809 

0.0827529 

SiO 

0.0275961 

0.00762716 

0 

0 

Na 

0.00589145 

0.00217786 

0.00599984 

0.00217104 

CO 

0.0055439 

0.262586 

0 

0 

SiCl 2 

0.00286999 

0.00102598 

0.00292221 

0.00100409 

SiCl 3 

4.49xl0“ 5 

3.47X10" 5 

4.61 xlO" 5 

4.41xl0" s 

SiCl 

7.08x10"® 

1.17X10" 6 

7.14X10" 6 

8.81 xlO' 7 

Cl 

5.91 xlO"® 

1.28xl0“ 5 

5.97X10" 6 

1.66xl0" 5 

Si 

2.43x10“® 

1.86xl0 -7 

2.43x10 r 6 

1.08xl0“ 7 

SiCl 4 

1.19x10“® 

1.99X10" 6 

1. 23x10“® 

3.28x10"® 

Naz 

4.44xl0~ 7 

6.07xl0 -8 

4.61xl0" 7 

6.04x10"* 

co 2 

6.20xl0“ 8 

1.06xl0 -5 

0 

0 

Si0 2 

8.86x10"® 

8.86x10"® 

0 

0 

Cl 2 

4.61 xlO" 10 

2.16x10"® 

4.69xlO" 10 

3.64x10"® 

NO 

1.24x10"“ 

3.87x10"“ 

0 

0 

NaO 

2.13xl0" 12 

2.85xl0" 12 

0 

0 

SiC 

2-lOxlO" 13 

2.10xl0 -13 

2.10xl0" 13 

2.10xl0" 13 

N 2 0 

5.18xl0" 16 

1.40xl0" 15 

0 

0 

°2 

2.95xl0' 18 

3.86x10"“ 

0 

0 

no 2 

1.15xl0“ 22 

1.31X10 -21 

0 

0 

Condensed 





SiC 

1 

1 

1 

1 

Si 

1 

0.0442988 

1 

0.0442988 

Si0 2 

1 

1 

0 

0 

NaCl 

0.570859 

0.456512 

0.586626 

0.591188 

Na 2 Si 2 O s 

0.00117705 

0.00058196 

0 

0 

NajSiOj 

0.00014994 

7.41xl0" 5 

0 

0 

C 

0.0442988 

1 

0.0442988 

1 
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TABLE VI.— CALCULATED ACTIVITIES FOR GASEOUS AND 
CONDENSED SPECIES WITH DIFFERENT CONDENSED 
PHASES PRESENT AT 1600 °C (2912 °F) 


Species 



SiC + NaCl + Si 

SiC + NaCl + C 

n 2 

0.457521 

0.361458 

0.52122 

0.523599 

NaCl 

0.370224 

0.304023 

0.413905 

0.420039 

SiO 

0.0904799 

0.043057 

0 

0 

N^Clj 

0.0384404 

0.1259221 

0.0480464 

0.0494808 

CO 

0.0277116 

0.25715 

0 

0 1 

Na 

0.010414 

0.00559088 

0.0112201 

0.00459793 

SiCl 2 

0.00508322 

0.00269329 

0.00547338 

0.00215791 

SiClj 

6.73xl0" 5 

5.45xl0 -5 

7.52xlO” 5 

7.34x1 0” 5 

SiCl 

2.57xl0" 5 

8.91X10 -6 

2.67X10” 5 

4.25X10” 6 

Cl 

1.72xl0 -5 

2.62X10 -5 

1.78xl0” 5 

4.41 xlO -5 

Si 

1.22xl0" s 

2.76x1 0 -6 

1.22X10” 5 

7.82xl0” 7 

SiCl 4 

1.25X10- 6 

1.55X10 -6 

1.45X10" 6 

3.51X10” 6 

Na 2 

1.03X10" 6 

2.98xl0' 7 

1.20X10” 6 

2-OlxlO” 7 

co 2 

5.94xl0 -7 

1.16xl0 -5 

0 

0 

SiOj 

7.01xl0” 8 

7.01xl0 -8 

0 

0 

Cl 2 

1.55xl0” 9 

3.63xl0 -9 

1.67X10” 9 

1.02xl0” 8 

NO 

8.49x10”” 

1.59xlO -10 

0 

0 

NaO 

1.80x10”” 

2.03xl0 -11 

0 

0 

SiC 

3.46xl0” 12 

3.46X10” 12 

3.46xl0” 12 

3.46xl0” 12 

n 2 o 

3.29xl0” 15 

5.47xl0” 15 

0 

0 

°2 

8.19xl0” 17 

3.62xlO” 16 

0 

! 0 

no 2 

3.39xl0“ 21 

1.33xlO -20 

0 

0 

Condensed 





Si 

1 

0.0642883 

1 

0.0642883 

SiC 

1 

1 

1 

1 

Si0 2 

1 

1 

0 

0 

NaCl 

0.32423 

0.266253 

0.362484 

0.367855 

Na,Si,0< 

0.00109197 

0.00066137 

0 

0 

Na ? SiO, 

0.00014874 

9.01 xlO -5 

0 

0 

C 

0.0642883 

1 

0.0642883 

1 
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TABLE Vn.— CALCULATED ACTIVITIES FOR GASEOUS AND 
CONDENSED SPECIES WITH DIFFERENT CONDENSED 
PHASES PRESENT AT 1700 °C (3092 °F) 


Species 


SiC + NaCl + 
Si0 2 + C 

SiC + NaCl + Si 

SiC + NaCl + C 

n 2 

0.311134 

0.272915 

0.50861 

0.511961 

NaCl 

0.275056 

0.24449 

0.434888 

0.445956 

SiO 

0.260651 

0.203483 

0 

0 

CO 

0.1206224 

0.253526 

0 

0 

Na 

0.0144191 

0.0113123 

0.0195785 

0.0090755 

NajCL, 

0.0108282 

0.00855532 

0.0270688 

0.0284641 

SiCl 2 

0.00704177 

0.00550906 

0.00954812 

0.004297 

SiCl 

7.46xl0“ 5 

5.1Sxl0 -5 

8.68xl0“ 5 

1.77xl0“ 5 

SiCI 3 

7.32x1 0 -5 

6.49xl0" 5 

0.00011563 

0.00011512 

Si 

5.17X10" 5 

3.15xl0“ 5 

5.17xl0“ 5 

4.75x10“* 

Cl 

4.07X10 -5 

4.62x1 0* 5 

4.74xl0“ 5 

0.00010496 

co 2 

4.67xl0“* 

1.26xl0“ 5 

0 

0 

Na, 

1.52x10“* 

9.34x1 0“ 7 

2.80x10“* 

6.01 xlO" 7 

SiCl 4 

9.07xl0' 7 

9.11xl0' 7 

1.67X10 -6 

3.67x10“* 

S 1 O 2 

4.48x1 0 -7 

4.48xl0” 7 

0 

0 

CI 2 

3.83xl0“ 9 

4.91xl0" 9 

5.19xl0“ 9 

2.54.xl0“ 8 

NO 

4.20x10“’° 

5.04xl0" 10 

0 

0 

NaO 

9.95x10“" 

l.OOxlO -10 

0 

0 

SiC 

4.28x10“" 

4.28x10“" 

4.28x10“" 

4.28x10“" 

n 2 o 

1.33X10" 14 

1.50xl0“ 14 

0 

0 

0 2 

1.64xl0“ 15 

2.68xl0“ 15 

0 

0 

NO, 

6.22xlO“ 20 

9.56xlO -20 

0 

0 

Condensed 





Si 

1 

0.091961 

1 

0.091961 

SiC 

1 

1 

1 

1 

Si0 2 

1 

1 

0 

0 

NaCl 

0.143584 

0.127628 

0.227019 

0.232796 

Na2Si 2 0 5 

0.00068369 

0.00053903 

0 

0 

NaoSiO, 

0.00010049 

7.92x1 0“ 5 

0 

0 

C 

0.091961 

1 

0.091961 

1 
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TABLE VUI.— CALCULATED 
DIAMETERS FOR 
MOLECULAR SPECIES 


Species 

Calculated 

diameter, 

(A) 

At 

3.542 

Cl 

3.613 

o 2 

3.467 

N 2 

3.798 

CO 

3.690 

SiO 

3.374 

SiCl 

3.748 

SiCl 2 

a 4.49 1 

SiClj 

a 5.234 

SiCl 4 

5.977 


'‘Estimated by analogy to other 
species; linear interpolation. 


TABLE IX.— CALCULATED MEAN FREE PATHS FOR SiCl 2 IN N 2 


Partial 

pressure 


Calculated mean free path (pm) 


(atm) 

n 2 

1200 °C 

1300 °C 

1400 °C 

1500 °C 

1600 °C 

1700 °C 

0.01 

37.2 

T1 

42.3 

44.8 

mm 

a 

0.07 

5.33 

mm 

6.0 

6.4 

■B 

WEM 


TABLE X.— CALCULATED DIFFUSIVTTY FOR SiCl 2 IN N 2 


Total 

pressure 

Calculated diffusivity (cm 2 /s) 

(atm) 

1200 °C 

1300 °C 

1400 °C 

1500 °C 

1600 °C 

1700 °C 

0.01 

73.9 

81.6 

89.5 

97.6 

106.0 

114.6 

0.07 

10.6 

11.7 

12.8 

13.9 

15.1 

16.4 


TABLE XI. — CALCULATED PARABOLIC RATE CONSTANTS kp FOR PINHOLE 

GROWTH 


Total 

pressure 

(atm) 

Calculated parabolic constants kp (cm 2 /s) 

1200 °C 

1300 °C 

1400 °C 

1500 °C 

1600 °C 

1700 °C 

0.01 

0.07 

4.33xl0 -7 

6.19xl0 -8 

1.88X10 -6 

2.68xl0” 7 

6.89X10 -6 

9.85xl0' 7 

1.67xl0" 5 

2.39x10 

3.72xl0 -5 

5.31X10 -6 

7.59xl0 -5 

1.08X10 -5 
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Figure 1 . — Electron micrograph of a polished cross section of 
the oxidation protected RCC material. 


o 2 CO 



Figure 2. — Schematic for Si0 2 formation on SiC. 
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Weight change, mg /cm 2 Weight change, mg 


2.0 




Figure 3. — Parabolic growth curve for Si0 2 forming on SiC at 1300 K, 1 atm oxygen. The parabolic 
rate constant is extracted from the lower plot. 
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log Psi„O y , bar Va P or pressure, P, bar 



log Po 2 . bar 

Figure 4. — Stability diagram for Si-C-0 at 1500 K. 5 



Temperature, T, K 


Figure 5. — Total vapor pressure from CO(g) and SiO(g) for 
carbon saturated SiC and Si0 2 , silicon saturated SiC and 
Si0 2 , and SiC and Si0 2 forming SiO and CO in a 3:1 molar 
ratio. 5 



Figure 6. — Volatility diagram for the Si-O system at 1 800 K (ref 1 0). 



Figure 7. — Vapor species over Na 2 0.2 (SiO 2 ) in vacuum as 
calculated via a free energy minimization computer code. 



Figure 8.— Vapor species over Na 2 0.2 (Si0 2 ) in .01 atm oxygen 
as calculated via a free energy minimization computer code. 
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Hangwire 


Electrobalance 




(b) 

Figure 9. — Schematic of balance apparatus for measuring vaporization rates, 
(a) Into a vacuum, (b) Into a gas atmosphere. 
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Figure 10. — Photomicrograph of thermally shocked RCC material, 
showing cracks in coating. 



Log P 02 

Figure 1 1 . — Phase stability for SiC-C-O system at 
2000 K (3141 °F). 



Temperature,°C 


Figure 1 2. — Plot of vapor pressures for volatile species in 
equilibrium with Si0 2 + SiC + NaCI (passive reaction). 



Figure 1 3. — Plot of vapor pressures for volatile species in 
equilibrium with SiC + NaCI (active reaction). 
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SiCl 2 + 0 2 — ►Si0 2 + 2CI 



' — Dense Si0 2 
/ scale 


Porous Si0 2 
product 


SiC 

SiC + 2CI — SiCI 2 + C 


(a) 



C/C 

(b) 


C/C 


(c) 



Figure 14. — Schematic representation of chloridation/ 
oxidation reaction mechanism, (a) Contamination of surface 
cracks, (b) Transient passive reaction in a pinhole, (c) Steady 
-state active reaction with pinhole growth and fume 
deposition at external surface. 



Figure 1 5. — Model for formation of a pinhole in a SiC coating on a 
C/C composite, (a) Schematic of a pinhole in a SiC coating on a 
C/C composite, (b) Partial pressure gradient in a pinhole. 



Figure 1 6. — Kinetics of formation of a pinhole in a SiC coating 
on a C/C composite. 
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APPENDIX A— ACTIVE TO PASSIVE AND PASSIVE TO ACTIVE TRANSITIONS 


Wagner’s paper on the active to passive transition (ref. 1) for silicon is the classic work in the field and can be 
adapted to SiC as is summarized here. Let us begin with a bare SiC surface. At very low oxygen pressures, a passive silica 
film does not form on an initially clean surface. Any oxygen molecules which do strike the surface react to form SiO(g) 
and CO(g) as follows: 


SiC(s) + 0 2 (g) = SiO(g) + CO(g) (Al) 

The limitation in the arriving oxygen flux which permits and retains active reaction can arise from two limiting 
conditions: 1. A partial vacuum where (he oxygen pressure approximates the total pressure, or 2. In the presence of a high 
pressure of an inert diluent gas, such that the dilute oxygen component can only reach the surface via gaseous diffusion 
through a stagnant boundary layer. Clearly, as long as SiC remains bare, the higher the oxygen flux, the more rapid this 
"active" reaction, as shown in figure (Al). The flux of Si loss in figure (Al) is simply proportional to the rate of oxygen 
arrival and reaction according to equation (Al). Obviously SiC(s) can be rapidly consumed under these conditions. 

As the pressure of oxygen is increased, an oxygen arrival flux is reached at which a stable film of Si0 2 (s) will 
form, causing an abrupt drop in rate shown in figure (Al). In the literature (ref. 2) the condition to form a stable scale is 
the equilibrium between the three solid phases SiC, Si, and Si0 2 , according to this expression: 

2 SiC(s) + Si0 2 (s) = 3 Si(c) + 2 CO(g) (A2) 

To attain this equilibrium and a "passive" condition, we need a minimum CO(g) pressure, which can be calculated 
from the standard Gibbs free energy change for Reaction (A2). Thus once reaction (Al) generates enough CO pr essure to 
achieve equilibrium for Reaction (A2), a stable SiQ 2 scale is formed and the evaporation loss is virtual ly stifled. It is 
aggnmwl that there are no kinetic barriers to these reactions. In its original publication the model presumed that the oxygen 
arrives at the surface by diffusion through some inert gas of much higher pressure, e.g. Ar, N 2 , or He. The flux of oxygen 
molecules, Jq^, striking a bare SiC surface may be written as follows: 

Jo 2 = [Do 2 /So 2 ] (Po 2 g -Po 2 WT (A3) 

Here 8 is the boundary layer thickness, P s the gas pressure in the bulk gas (far side of the boundary layer), P 1 the pressure 
at the gas/SiC interface (generally assumed to be negligible), R the gas constant, and T the absolute temperature. The flux 
of CO(g) molecules, J co , leaving the surface is: 

J CO = [D C o/8 co ] (P C o - Pco^/RT < A4 > 

Now P co g is assumed to be negligible. The fluxes are shown schematically in figure (A2). An analogous expression to 
(A4) can be written for the SiO flux. However we need not consider this, since it does not enter into the stable Si0 2 /SiC 
equilibrium condition (equation (A2)). The fluxes can be equated according to equation (Al), Jq 2 = Jco = ^SiO 
result is: 

Po 2 = (D C0 /D 02 ) (5o 2 /8 co ) Pcq* (A5) 

This describes the steady state value for P C0 ‘ for active oxidation. As Wagner pointed out, the ratio of boundary layer 

thicknesses can be approximated as (Dq /D^q) 1 ^. Thus the necessary Pq for the active-to-passive transition becomes: 

2 

Po 2 >(D co /D 02 ) 1/2 Pco 6 " (A6) 

where P co a * is the equilibrium CO pressure calculated for Reaction (A2). The ratio, D co /D 0 (as inferred from the 
molecular diameters of table VIII in the mam text) is close to one. A plot of the predicted acfive-to-passive transition 
temperatures is shown in figure (A3). Measured values are within an order of magnitude of this and indicate the proper 
temperature dependence (ref. 2). 

Suppose, however, that one starts with a passive Si0 2 (s) scale and the oxygen pressure is gradually reduced. At 
what point does the scale become unstable and active oxidation occur? Clearly, the above criterion does not hold and we 
must invoke a new criterion, which leads to the observed hysteresis, as indicated in figure (Al). 
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The factors which lead to the breakdown of passive Si0 2 (s) need to be identified. Note that for this analysis, any 
consideration of a NaCl deposit is avoided. We have an oxygen flux striking the surface, which contributes to scale growth. 
We also have SiO(g), Si0 2 (g), and 0 2 (g) fluxes involved in the following vaporization processes: 

Si0 2 (s) = SiO(g) + 1/2 0 2 (g) (A7) 

Si0 2 (s) = Si0 2 (g) (A8) 

At low oxygen partial pressures, reaction (A7) dominates (ref. 3) and we can ignore the flux of Si0 2 (g). The flux of 0 2 (g) 
is given by: 


Jo 2 =[Do 2 / 8 o 2 ](P o| - p o 2 ^ T (A9) 

Here the terms have the same definitions as given previously, with Pq 2 being the pressure of 0 2 at the surface. In this 

instance, neither P 0 g or P 0 S can be neglected. The flux of SiO(g) leaving the surface is given as: 

2 2 

J SiO = t D Sic/8siol ^Sic/ * P SiO S ) /RT (A10) 

The two fluxes are shown schematically in figure (A4) and are related by the expression J Q = 1/2 J Si0 . Now we can 
neglect P Si0 g and express P SiQ s in terms of P 02 s via the equilibrium constant for Reaction (A7): 

K = P Si o S ( p o 2 S ) 1/2 (A11) 

Combining eq uati ons (A9) to (All) and correlating the boundary layer thickness to diffusivities, as done previously, leads 
to: 

2(P 0 g 2 - P 0 S 2 > = K (D Si( yDo 2 ) 1/2 /(Po S 2 ) m (A12) 

As explained by Wagner (ref. 1), the left-hand side of equation (A12) is proportional to the supply of oxygen and the right- 
hand side is proportional to the consumption of oxygen. A stable Si0 2 scale is maintained only when the oxygen supplied 
is equal to or gr eate r than the amount consumed. The minimum P 0 g for which the conditions of equation (A12) is met can 
be calculated by expressing equation (A12) as a function of P 0 g and calculating the m i n i mu m for that expression: 

f(P Q g ) = 2P 0 s 2 + K (D Si( yD 02 ) 1/2 /(P 0 s 2 ) m = 2 P Q g (A13) 

The minimum value is obtained by taking the derivative of equation (A13) with respect to P 0 S 2 and setting the resultant 
expression equal to zero: 

P 0 S 2 = (KM) 2 * (D Si0 /D 02 ) 1/3 (A14) 

The value of P 0 g can be calculated by combining equations (A14) and (A12): 

P Q g = 1/2 [2/(4)2* + (4) 1 *] K 2 * (D Si0 /D 02 ) 1/3 (A15) 

This is the oxygen pressure at which active oxidation begins; such values have been plotted onto figure (A3) which shows 
the results of this calculation for a range of temperatures. Note the large difference between the active-to-passive and the 
passive-to-active transition points. Also, recall the outer layer on the RCC panels is likely to be a silica-satured sodium 
silicate, not pure silica. So these boundaries are only approximations. 

Since the RCC panels start with a coating of Si0 2 (s), it is this lower transition pressure that induces active behavior 
that concerns us. Clearly the transition to true active oxidation is only expected at this extreme condition of very low 
pressure. 
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Figure A1 . — Schematic showing active to passive and passive to 
active transitions as a function of PfOg). 



1/T(K) 


Figure A3. — Calculated active-to-passive and passive-to active 
transitions for SiC as a function of temperature. 




Boundary Layer 

Figure A2. — Schematic showing diffusion of O 2 , 
SiO, and CO through the boundary layer at a bare 
SiC surface. Boundary layers for each gas are 
shown as the same thickness due to close 
molecular weights. 



SiOJ 



Boundary Layer 


°2 


SiO 


Figure A4. — Schematic showing diffusion of 
O 2 through the boundary layer away from 
Si02- Boundary layers for each gas are shown 
as the same thickness due to close molecular 
weights. 
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APPENDIX B— OXIDATION OF C/C BENEATH A SiC COATING CONTAINING PINHOLES 


When pinholes form in a SiC coating, the underlying C/C composite substrate is exposed to an oxidizing 
environment through the pinholes, and oxidation of the carbon occurs. A calculation of the oxidation rate though a flawed 
SiC coating on the C/C composite is made in the following. Some parameters used are as follows: 

(1) Thickness of the SiC coating: 20 to 40 mils (508 to 1016 pm) 

(2) Diameter of the pinhole: 20 to 40 mils (508 to 1016 pm) 

(3) Environment: 001 to 0.07 atm-air (pressure created by shock wave) 

(4) Temperature: 1473 to 1973 K 

The oxidation of C layers in composites has been reported in references 1 and 2, and a similar treatment for the 
oxidation of TaC is treated in reference 3. A modification of the model in reference 1 is employed for oxidation of a C/C 
composite underneath a pinhole in a SiC coating, as shown in figure (Bl). The pinhole is filled with four gaseous species: 
inert N 2 , the active oxidant 0 2 from air, and the two oxidation products CO and C0 2 . At the C/gas phase interface, the 
oxygen partial pressure is very low and the primary C oxidation product at very high temperature is CO. Under these 
conditions CO and 0 2 cannot coexist. Hence we have two regions in the pore— one with 0 2 and C0 2 and one with CO and 
C0 2 . The CO forms by reaction of C with C0 2 at the base of coating, diffuses outward and at a more forward site in the 
pinhole, Xf, it is oxidized a gain to form C0 2 . The C0 2 oxidation product diffuses both inward as an oxidant to oxidize C at 
x = 0 and outward as a product through the pinhole to the environment. The gas within the voids in C/C caused by 
oxidation through the pinhole is assumed to remain in equilibrium with pure C (virtually pure CO). The concentration 
(partial pressure) profiles of the gaseous species in the pinhole are shown in figure (B2). The two regions (I and II) are 
identified in the pinhole, as labelled in figure (B2). The concentration profiles and the rate of carbon oxidation can be 
obtained through analysis of the diffusion fluxes for the gaseous species in these two regions. 

(I) Diffusivities of 0 2 , CO and C0 2 in N 2 

To calculate the oxidation rate of carbon at the base of a pinhole, diffusivities for 0 2 , CO and C0 2 in N 2 are 
required. These diffusivities are calculated by the following equation (ref. 4): 

D ; . N2 = (2/3)(k B /re) 3/2 [l/(2m i ) + l/(2m N2 )] I/2 T 3/2 /P T /[(d 1 + d N2 )/2] 2 (Bl) 

where k B is the Boltzmann constant, ni; the weight of species i per molecule, P T the total pressure, T the absolute 
temperature, and dj the diameter of species i. The calculated diffusivities are listed in table Bl. 


TABLE Bl— DIFFUSIVITIES (cm 2 /s) OF O z , CO AND 
C0 2 IN N 2 * 


Species, 

dj(A) 5 

o 2 , 

3.467 

CO, 

3.690 

C0 2 

3.941 

P(atm) 

0.01 

0.07 

0.01 

0.07 

0.01 

0.07 

1473K 

116.4 

16.6 

113.2 

16.2 

95.8 

13.7 

1S73K 

128.5 

18.4 

124.9 

17.8 

105.8 

15.1 

1673K 

140.9 

20.1 

137.0 

19.6 

116.0 

16.6 

1773K 

153.7 

22.0 

149.4 

21.3 

126.6 

18.1 

1873K 

166.9 

23.8 

162.3 

23.2 

137.4 

19.6 

1973K 

180.5 

25.8 

175.4 

25.1 

149.6 

21.2 


♦Diameter of N 2 is 3.798 A (ref. 5). 


(II) Mass Transport in Region I 

In region I of figure B2, a steady-state mass transport calculation is provided as follows. The inward flux of 0 2 is 
independent of distance in region I and given by: 
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J'i= -DjCdCj/dx) + vC ; 


(B2) 


where D ; and Cj are the diffusivity and concentration of species i, and v is the local molar average velocity, which is 
defined as: 

v = KC^/Sq = Kty/Cp (B3) 


where Cj, v ; , and are the concentration, local molar velocity and the flux of species i, respectively, and C T is the total 
concentration of gaseous species. The net flux of N 2 is zero in both the regions I and II because neither a sink nor a source 
for N 2 exists in the pinhole. The concentration and flux of CO are negligible in region I. From the mass balance, the 
inward flux of 0 2 equals the outward flux of C0 2 in region I, i.e.. 


JV = - J r 


CO, 


(B4) 


and by inserting equation (B4) into equation (B3), v is equal to zero. The concentration profiles for both 0 2 and C0 2 are 
linear because of the constant dCj/dx in equation (B2). From equation (B2) for 0 2 and C0 2 : 


Jo 2 = -Do 2 C 0 ° 2 /(L-x f ) 


(B5) 


and 

Jco 2 = ^C0 2 ^£o 2 ^ “ 

From equations (B4) and (B6), 

C C0 2 = C 0 2 D 02 /D C o 2 <® 7 > 

where Cq is the oxygen concentration at the external surface, and CJ^. is the C0 2 concentration at x = Xf. The 
concentration profiles for 0 2 and C0 2 in region I are plotted in figure (B2). 


(Ill) Mass Transport in Region 13 

From a mass balance, the following relations are obtained: 

j N 2 = 0 

Jco = -2 Jco 2 

and inserting equations (B8) and (B9) into equation (3) yields: 

v = -Jco 2 /C t 

The flux of C0 2 in the region II is given by: 

Jco 2 = -Dco,^/^) “ ( -'C0 2 ^C0 2 ^T 


(B8) 

(B9) 

(BIO) 


(Bll) 


Rearrangement of equation (Bll) gives: 


•^CO^ ~ -[C t Dco 2 /(Ct + C C o 2 )]dC C o 2 (B12) 

Equation (B12) is integrated from x = 0 to x^ on the left-hand side and from Cco 2 = ^*C0 2 ^ C0 2 on right-hand 

side, where the superscript double star denotes the reaction side at x = 0. The result is: 

Jco 2 = -(C T D c02 /x f )ln[(C x + C £o2 )/(C x + C*^)] (B13) 
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Because C*^ 0 is much lower than C£ 0 . equation (B13) is simplified as: 

2 2 

Jco 2 = -(CtDco^MI + C* CQ2 /C t ] = -(OPco^) ln[l + CS 2 D 02 /(C t D C02 >] (B14) 

A mass balance between CO and 0 2 results in the following relation: 

Sqo = -2Jq 2 (B15) 

Combining equations (B9) and (B15) yields: 

J 0 2 = J co 2 

From equations (B5), (B14) and (B16), 

-(CTDco^lntl + Cg 02 D 02 /(C T D c02 )] = - C^D 02 /(L - xj (B17) 

Equation (B17) can be rearranged to give: 

L/Xf = 1 + CS 2 D 02 /(C T D c02 )/ln[l + ^Do^D^)] (B18) 


According to equation (Bll), D^/D cc>2 is independent of temperature. For a system with a given C^Cp, (i.e., Pq^j) 

the environment, the ratio of the lengths of the regions I to II is a constant For the environment of air, Pq^t is °- 21 311(1 

the calculated value for L/Xf is 2.123. 

Similarly, the flux of CO in the region II can be obtained as follows: 

Jco = -D co (dC co /dx) + C co (J£o + Jco 2 )/C x (B19) 

From equation (B9), equation (B19) can be rewritten as: 

J£ 0 dx = -[2 C x D co /(2C t - C co )]dC co (B20) 

Integration of equation (B20) from x = 0 to x f and C co = C*^ 0 to C£ 0 yields: 

Jco = -(2C T D co /x f )ln[(2C T - C* c * 0 )/(2C T - C&>)] (B21) 

Because C£ 0 is much less than C*<* 0 , equation (B21) is simplified as: 

J£ 0 = -(2C T D co /x f )In[ 1 - C* c ? 0 )/(2C t )] , (B22) 

By combining equations (B9), (B14) and (B22), C^q is obtained as: 

C* c *o = 2C T { 1 - [1 + CS 2 D 02 /(C x D C o 2 )]- d C02®C0} (B23) 

The concentration profiles of both CO and C0 2 are nonlinear in the region II, as shown in figure (B2). 

(IV) Oxidation Rate for C/C 

The weight-loss rate due to C/C oxidation per unit pinhole area in a SiC coating is related to the flux of C0 2 , according to 
a mass balance, as follows: 

dW c /dt = U c r C o 2 = M c Jq 2 = M c CS 2 Do 2 /(L - x*) = M c P 02 D 02 /[(L - Xf)RT] (B24) 

The calculated weight-losses per cm 2 of pinhole area are the same for P = 0.01 and 0.07 atm, as plotted in figure (B3). The 
weight-loss rate per pinhole is given by: 
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(B25) 


W c = ^(dWc/dt) = i^tcMcJ^ = r 2 nM c C^ 2 D 02 /(L - Xf) 

where r is the radius of the pinhole. The calculated results for pinhole of diameters of 20 to 40 mils at temperatures of 
1473 and 1973 K are shown in figures (B4) and (B9). 

If a hemispherical cavity is assumed to grow at the bottom of a pinhole, the weight loss rate of C can be converted 
to the size of a cavity as follows: 


W c = dc(dV/dt) = dc[d(2iti^/3)dt] (B26) 

where V is the volume of the cavity, d c the density of C/C, and r c the radius of the cavity. Integration of equation (B26) 
yields: 

r c = [3W c t/(27cdc)] 1/3 = {3r 2 M c P 02 D 02 t/[RT(L - x f )(2d c )] } 1/3 (B27) 

The density of C is in a range of 1.8 to 2.25 g/cm 3 (ref. 6), and 2 g/cm 3 is chosen in the following calculation. The 
calculated results are plotted in figures (BIO) and (Bll) for temperatures of 1473 to 1973 K and different pinhole sizes. 

The size of a cavity formed is strongly dependent on the size of a pinhole, and weakly dependent on temperature, as shown 
in figure (BIO). 

(V) Discussion 

(i) Ratio of L/Xf 

The ratio of L/Xf is very important to the oxidation rate of C because the oxidation rate is proportional to the flux of 
0 2 , which is inversely proportational to L - Xf. As analyzed in the section III, the ratio of L/Xf is independent of both 
temperature and total pressure if P Q /P T is fixed. But the ratio Pq^t plays an important role in the oxidation rate by 
changing the ratio of L/x^ At a given temperature and a given P 0 , the ratio of L /Xf decreases with increasing P T , 
according to equation (B18). Two extreme cases are discussed in me following: 

(a) P T » Pq 2 

In this case, L/x^ = 2 because ln[l + C^D^/fCpDco^] is approximately equal to ^ 2 D 0 2 ^ C T D C 0 2 ) in equation (B18). 
Also, minimum value of L/Xf is 2. 

(b) Pj = Pq 2 

This is a case of oxidation in a pure 0 2 environment Equation (B18) becomes: 

L/Xf = (1 + D 02 D C02 )/ln(l + D 02 /D c02 ) = 1 + l/ln2 = 2.44 (B28) 

if D 0 (in C0 2 ) is assumed equal to D co (in CO), which is the maximum value of L/Xf. 

2 2 

(ii) Effect of Total Pressure on Oxidation Rate 

The diffusivity of gaseous species is inversely proportional to P T , as shown in equation (Bl). For a system with a 
given Pq 2 at a given temperature, the oxidation rate is inversely proportional to P T , the same as Dq 2 _ 

Because Cq 2 is proportional Pq 2 and Dq 2 is inversely proportional to P T , the product of C 02 D 02 is unchanged if 
Pq 2 /P t is fixed. For a given Pq 2 /P t and a given temperature, the oxidation rate is unchanged even if both Pq 2 and P T are 
changed. 
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(iii) Temperature Effect on Oxidation Rate 


According to the above calculation, the ratio of L/Xf is independent of temperature. For a system with a given 
Po 2 /P t , the oxidation rate is proportional to Dq 2 , which is proportional T 1 ' 5 . In other words, the oxidation rate for a given 
pinhole population increases with the square root of temperature if the ratio of Pq^/P -j- is fixed. The growth rate of a cavity 
in C/C at the bottom of a pinhole is proportional to T 1/6 , according to equation (B27), and as temperature increases from 
1473 to 1973 K, the growth rate only increases 5 percent The growth rate is virtually dependent on the size of a cavity if 
the ratio of Pq 2 /P t is fixed. 
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